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ABSTRACT 

Context. The mechanisms giving rise to diffuse radio emission in galaxy clusters, and in particular their connection with cluster 
mergers, are still debated. 

Aims. We aim to obtain new insights into the internal dynamics of the cluster Abell 1240, showing the presence of two roughly 
symmetric radio relics, separated by ~ 2 fcjjj Mpc. 

Methods. Our analysis is mainly based on redshift data for 145 galaxies mostly acquired at the Telescopio Nazionale Galileo and on 
new photometric data acquired at the Isaac Newton Telescope. We also use X-ray data from the Chandra archive and photometric 
data from the Sloan Digital Sky Survey (Data Release 7). We combine galaxy velocities and positions to select 89 cluster galaxies and 
analyze the internal dynamics of the Abell 1237 + Abell 1240 cluster complex, being Abell 1237 a close companion of Abell 1240 
towards the southern direction. 

Results. We estimate similar redshifts for Abell 1237 and Abell 1240, <z) = 0. 1935 and (z) = 0. 1948, respectively. For Abell 1237 we 
estimate a line-of-sight (LOS) velocity dispersion tr v ~ 740 km s -1 and a mass M ~ 6 xlO 14 thjg M Q . For Abell 1240 we estimate a 
LOS <r v ~ 870 km s -1 and a mass range M ~ 0.9 —1.9 xlO 15 fci M Q , which takes into account its complex dynamics. Abell 1240 is 
shown to have a bimodal structure with two galaxy clumps roughly defining the N-S direction, the same one defined by the elongation 
of its X-ray surface brightness and by the axis of symmetry of the relics. The two brightest galaxies of Abell 1240, associated to the 
northern and southern clumps, are separated by a LOS rest-frame velocity difference Vrt ~ 400 km s -1 and by a projected distance 
D ~ 1.2 /!,„' Mpc. The two-body model agrees with the hypothesis that we are looking at a cluster merger occurred largely in the 
plane of the sky, with the two galaxy clumps separated by a rest-frame velocity difference Vrt ~ 2000 km s at a time of 0.3 Gyrs 
after the crossing core, while Abell 1237 is still infalling onto Abell 1240. Chandra archive data confirm the complex structure of 
Abell 1240 and allow us to estimate a global X-ray temperature T x = 6.0 + 0.5 keV. 

Conclusions. In agreement with the findings from radio data, our results for Abell 1240 strongly support the "outgoing merger shocks" 
model to explain the presence of the relics. 

Key words. Galaxies: clusters: individual: Abell 1240, Abell 1237 - Galaxies: clusters: general - Galaxies: kinematics and dynamics 



1. Introduction 

Merging processes constitute an essential ingredient of the evo- 
lution of galaxy clusters (see Feretti et al. 2002b for a review). 
An interesting aspect of these phenomena is the possible connec- 
tion of cluster mergers with the presence of extended, diffuse ra- 
dio sources: halos and relics. The synchrotron radio emission of 
these sources demonstrates the existence of large-scale cluster 
magnetic fields and of widespread relativistic particles. Cluster 
mergers have been suggested to provide the large amount of en- 
ergy necessary for electron reacceleration up to relativistic en- 
ergies and for magnetic field amplification (Feretti 1999; Feretti 
2002a; Sarazin 2002). Radio relics ("radio gischts" as referred 
by Kempner et al. 2003 ), which are polarized and elongated ra- 
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dio sources located in the cluster peripheral regions, seem to be 
directly associated with merger shocks (e.g., Ensslin et al.[1998; 
R oettige r et al. [19991 Ensslin & Gopal-Krishna [20TJT] Hoeft et 
al. 120041) . Radio halos, unpolarized sources which permeate the 
cluster volume similarly to the X-ray emitting gas, are more 
likely to be associated with the turbulence following a cluster 
merger (Cassano & Brunetti 2005). However, the precise radio 
halos/relics formation scenario is still debated since the diffuse 
radio sources are quite uncommon and only recently one can 
study these phenomena on the basis of a sufficient statistics (few 
dozen clusters up to z ~ 0.3, e.g., Giovannini et al. |1999t see also 
Giovannini & Feretti 2002; Feretti 2005a) and attempt a classi- 
fication (e.g., Kempner et al. 2003; Ferrari et al. 2008 ). 

There is growing evidence of the connection between dif- 
fuse radio emission and cluster merging, since up to now diffuse 
radio sources have been detected only in merging systems. In 
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most of the cases the cluster dynamical state has been derived 
from X-ray observations (see Buote 2002 ; Feretti 2006 and 2008 
and refs. therein). Optical data are a powerful way to investigate 
the presence and the dynamics of cluster mergers (e.g., Girardi 
& Biviano |20021 >. too. The spatial and kinematical analysis of 
member galaxies allow us to detect and measure the amount 
of substructure, to identify and analyze possible pre-merging 
clumps or merger remnants. This optical information is really 
complementary to X-ray information since galaxies and intra- 
cluster medium react on different time scales during a merger 
(see, e.g., numerical simulations by Roettiger et al. I1997I ). In 
this context we are conducting an intensive observational and 
data analysis program to study the internal dynamics of clusters 
with diffuse radio emission by using member galaxies (Girardi 
et al. l20"07l and refs. therein Q). 

During our observational program we have conducted an in- 
tensive study of the cluster Abell 1240 (hereafter A 1240). 

A 1240 is a very rich, X-ray luminous, Abell cluster: 
Abell richness class = 2 (Abell et al. [I9"g9l ; L x (0. 1-2.0 
keV)= 8.3xl 43 h~l erg s" 1 recovered from ROSAT data (David 
et al. 119991 correcting for our cluster redshift, see below). 
Optically, the cluster center is not dominated by any single 
galaxy - it is classified as Bautz-Morgan class III (Abell et al. 
[19891 

Kempner & Sarazin (2001) revealed the presence of two 
roughly symmetric radio relics from the Westerbork Northern 
Sky Survey. They appear to either side of the cluster center, north 
and south, at distances of ~ 6' and 7'. Kempner & Sarazin also 
noticed that A 1240 shows an elongated X-ray morphology (re- 
covered from ROSAT observations) consistent with a slightly 
asymmetric merger with the apparent axis roughly aligned with 
the axis of symmetry of the relics (see also Bonafede et al. 2009 ). 
The presence of double relics was confirmed by recent, deep 
VLA observations (Bonafede et al. 2009; see Fig. [T}. Very few 
other clusters with double relics have been observed: Abell 3667 
(Rottgering et al. \lWf\, Abel l 3376 (Bagchi et al. [2 006ll, Abell 
2345 (Giovannini et al. [19991 Bona fede et al. 120091 and RXCJ 
1314.4-2515 (Feretti et al. I2005bl Venturi et al. l2007l l. The 
relics of Abell 3667 were explained with the "outgoing merger 
shocks" model (Roettiger et al. [1999). Observations of Abell 
3376 agree with both the "outgoing merger shocks" and the "ac- 
cretion shocks" models (Bagchi et al. 2006). In the case of Abell 
2345 the observations are difficult to reconcile with theoretical 
scenarios (Bonafede et al. 2009). Instead, more data are needed 
for RX CJ 1314.4-2515 (Feretti et al. |2005bJ and Venturi et al. 
2007). As for A 1240, the detailed analysis of its radio properties 
is in agreement with the "outgoing merger shocks" (Bonafede et 
al. 120091 1, but the main global properties are unknown and the 
internal cluster dynamics was never studied. 

Indeed, few spectroscopic data have been reported in the 
field of A 1240 (see NED) and the value usually quoted in the lit- 
erature for the cluster redshift (z = 0.159; see, e.g., David et al. 
1999) is given by Ebeling et al. dl996l l, on the basis of the 10th- 
ranked cluster galaxy. The real cluster redshift, as estimated in 
this paper, is rather z = 0.195. Even poorer information is known 
for Abell 1237 (hereafter A 1237), a close southern companion 
of A1240, having richness class = 1 and Bautz-Morgan class III 
(Abell etal. [19891 . 

Recently, we performed spectroscopic and photometric ob- 
servations of the A1237+A1240 complex with the Telescopio 
Nazionale Galileo (TNG) and the Isaac Newton Telescope 
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(INT), respectively. Our present analysis is mainly based on our 
new optical data and X-ray Chandra archival data. We also use 
the few public redshifts and photometric data from the Sloan 
Digital Sky Survey (SDSS, Data Release 7). This paper is orga- 
nized as follows. We present our new optical data and the clus- 
ter catalog in Sect. 2. We present our results about the cluster 
structure based on optical and X-ray data in Sects. 3 and 4, re- 
spectively. Finally, we briefly discuss our results and give our 
conclusions in Sects. 5 and 6. 

Unless otherwise stated, we give errors at the 68% confi- 
dence level (hereafter c.l.). Results with a c.l. below 90% are 
considered very poorly/no significant. The values of these c.ls. 
are generally not explicitly listed throughout the paper. 

Throughout this paper, we use Ho = 70 km s Mpc in a 
flat cosmology with Qo - 0.3 and £2 a = 0.7. In the adopted cos- 
mology, 1' corresponds to ~ 194 h^J kpc at the cluster redshift. 

2. New data and galaxy catalog 

Multi-object spectroscopic observations of A 1240 were carried 
out at the TNG telescope in December 2006 and December 
2007. We used DOLORES/MOS with the LR-B Grism 1, yield- 
ing a dispersion of 187 A/mm. In December 2006 we used the 
old Loral CCD, with a pixel size of 15 fim, while in December 
2007 we used the new E2V CCD, with a pixel size of 13.5 
[im. Both the CCDs are matrices of 2048 x 2048 pixels. In to- 
tal we observed four MOS masks (2 in 2006 and 2 in 2007) 
for a total of 142 slits. We acquired three exposures of 1800 
s for each mask. Wavelength calibration was performed using 
Helium-Argon lamps. Reduction of spectroscopic data was car- 
ried out with the IRAF0 package. 

Radial velocities were determined using the cross- 
correlation technique (Tonry & Davis 1979) implemented in the 
RVSAO package (developed at the Smithsonian Astrophysical 
Observatory Telescope Data Center). Each spectrum was corre- 
lated against six templates for a variety of galaxy spectral types: 
E, SO, Sa, Sb, Sc, Ir (Kennicutt |1992| l. The template producing 
the highest value of %, i.e., the parameter given by RVSAO and 
related to the signal-to-noise ratio of the correlation peak, was 
chosen. Moreover, all spectra and their best correlation functions 
were examined visually to verify the redshift determination. In 4 
cases (IDs 55,61, 69, 129, 137;seeTable[B we took the EMS AO 
redshift as a reliable estimate of the redshift. Our spectroscopic 
survey in the field of A 1240 consists of spectra for 118 galaxies 
with a median nominal error on cz of 50 km s _1 . The nominal 
errors as given by the cross-correlation are known to be smaller 
than the true errors (e.g., Malumuth et al. 1992; Bardelli et al. 
[T9941 Ellingson & Yee [T994l Quintana et al. 2000; Boschin et al. 
2004). Double redshift determinations for four galaxies allowed 
us to estimate real intrinsic errors. We compared the first and sec- 
ond determinations computing the mean and the rms of the vari- 
able (zi - Zi)/ Jerrf + ertf- We obtained mean = 0.3 ±1.4 and 
rms = 2.8, to be compared with the expected values of and 1. 
The resulting mean shows that the two sets of measurements are 
consistent with having the same velocity zero-point. According 
to the ^ 2 -test the high value of the rms suggests that the errors 
are underestimated. Only when nominal errors are multiplied by 
a ~ 2 factor the rms is in acceptable agreement with the value 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. INT /?-band image of the cluster A1240 (North at the top and East to the left) with, superimposed, the contour levels of 
the Chandra archival image ID 4961 (thin contours; photons in the energy range 0.5-2 keV) and the contour levels of a VLA radio 
image at 1.4 GHz (thick contours; Bonafede et al 2009 ). Arrows show the positions of the two radio relics. Boxes highlight the 
brightest galaxies of A1240: BCGN, BCGS and BCGE (see text). 



of 1 . Therefore, hereafter we assume that true errors are larger 
than nominal cross-correlation errors by a factor 2. For the four 
galaxies with two redshift estimates we used the weighted mean 
of the two measurements and the corresponding errors. 

We also use 32 public redshift data as taken from NED 
within a box of 45x45' from the cluster center. They come from 
the SDSS (Data Release 7). Before to proceed with the merging 



between our and published catalogs we payed particular atten- 
tion to their compatibility. Five galaxies are in common between 
SDSS and TNG data. For them we computed the mean and 

the rms of the variable (z our - ziit)/ Jerrl ur + err^. We obtained 

mean - -0.1 + 0.6 and rms = 1.5, to be compared with the ex- 
pected values of and 1 . The resulting mean shows that the two 
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Fig. 2. INT /?-band image of the A1237+A1240 complex (North at the top and East to the left). Circles and annuli indicate cluster 
members of A1240 and A1237, respectively (see Table[TJ. Boxes indicate non-member galaxies. 



sets of measurements are consistent with having the same veloc- 
ity zero-point, and the value of rms is compatible with a value of 
1 according to the ^- 2 -test. Thus we added the redshifts coming 
from the literature. For the five galaxies in common we used the 
weighted mean of the two redshift determinations and the cor- 
responding errors. We obtained a final catalog of 145 galaxies 
with measured radial velocities. 

As far as photometry is concerned, our observations were 
carried out with the Wide Field Camera (WFC), mounted at 
the prime focus of the 2.5m INT telescope. We observed the 
A1237+A1240 complex with Bu and Ru in photometric con- 
ditions. The R image was obtained in December 19th 2004 with 



a seeing condition of 3.0". We got the B image in May 14th 2006 
with a seeing of about 1.1". 

The WFC consists of a four-CCD mosaic covering a 
33'x33' field of view, with only a 20% marginally vignetted 
area. We took nine exposures of 720 s in Z?h and 300 s in Ru 
Harris filters (a total of 6480 s and 2700 s in each band) devel- 
oping a dithering pattern of nine positions. This observing mode 
allowed us to build a "supersky" frame that was used to cor- 
rect our images for fringing patterns (Gullixson 1992). In addi- 
tion, the dithering helped us to clean cosmic rays and avoid gaps 
between the CCDs in the final images. The complete reduction 
process (including flat fielding, bias subtraction and bad-column 
elimination) yielded a final coadded image where the variation 
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of the sky was lower than 0.4% in the whole frame. Another 
effect associated with the wide field frames is the distortion of 
the field. In order to match the photometry of several filters, a 
good astrometric solution is needed to take into account these 
distortions. Using the imcoords IRAF tasks and taking as a ref- 
erence the USNO B 1 .0 catalog, we were able to find an accurate 
astrometric solution (rms ~ 0.5") across the full frame. The pho- 
tometric calibration was performed using Landolt standard fields 
obtained during the observation. 

We finally identified galaxies in our Bu and Ru images and 
measured their magnitudes with the SExtractor package (Bertin 
& Arnouts 1996 ) and AUTOMAG procedure. In a few cases (e.g. 
close companion galaxies, galaxies close to defects of the CCD) 
the standard SExtractor photometric procedure failed. In these 
cases we computed magnitudes by hand. This method consists 
in assuming a galaxy profile of a typical elliptical galaxy and 
scaling it to the maximum observed value. The integration of 
this profile gives us an estimate of the magnitude. This method 
is similar to PSF photometry, but assumes a galaxy profile, more 
appropriate in this case. 

We transformed all magnitudes into the Johnson- 
Cousins system (Johnson & Morgan 119531 Cousins 
[19761 . We used B — Bu + 0.13 and R — Rn 
as derived from the Harris filter characterization 
(http://www.ast.cam.ac.uk/~wfcsur/technical/photom/colours/) 
and assuming a B - V ~ 1.0 for E-type galaxies (Poggianti 
119971 1. As a final step, we estimated and corrected the galactic 
extinction Ab ~ 0.10, Ar ~ 0.06 from Burstein & Heiles's 
( 1982) reddening maps. 

We estimated that our photometric sample is complete down 
to R = 20.5 (22.0) and B = 22.0 (23.0) for S/N = 5 (3) within 
the observed field. 

We assigned magnitudes to all galaxies of our spectroscopic 
catalog. We measured redshifts for galaxies down to magnitude 
R ~ 20, but a high level of completeness is reached only for 
galaxies with magnitude/? < 19 (~45% completeness). 

Table [T] lists the velocity catalog (see also Fig. [2]): identi- 
fication number of each galaxy, ID (Col. 1); right ascension 
and declination, a and 5 (J2000, Col. 2); B and R magni- 
tudes (Cols. 3 and 4); heliocentric radial velocities, v = czq 
(Col. 5) with errors, Av (Col. 6); redshift source (Col. 7; T:TNG, 
S:SDSS); member assignment (Col. 8; 1:A1240, 2:A1237, 
: b ackground/f oreground) . 



Table 1. Velocity catalog of 145 spectroscopically measured 
galaxies in the field of the A1237+A1240 complex. In Col. 1, 
IDs 75, 56, 111, 29 and 27 (in boldface) highlight BCGS, 
BCGN, BCGE, BCG1 and BCG2, respectively (see text). 
Asterisks in Col. 8 highlight emission line galaxies. 



ID 




a,c5(J2000) 


B 


R 


V 

(km 


Av 


So. 


CI. 


1 


1 1 


22 00.12, 


+42 54 25.2 


17.76 


15.75 


23230 


53 


s 





2 


1 1 


22 03.36, 


+43 11 41.6 


18.53 


18.25 


180587 


338 


s 





3 


1 1 


22 04.44, 


+43 01 35.4 


18.57 


16.87 


37090 


40 


s 


o 


4 


1 1 


22 24.24, 


+42 53 08.9 


18.66 


17.46 


44576 


27 


s 





5 


1 1 


22 41.52, 


+43 00 41.4 


18.47 


17.17 


21546 


49 


s 


o 


6 


1 1 


22 46.20, 


+42 57 37.4 


19.28 


17.28 


33915 


46 


s 





7 


1 1 


22 59.88, 


+43 06 43.6 


18.48 


17.10 


34993 


32 


s 


o 


8 


1 1 


23 08.52, 


+43 10 41.9 


18.93 


17.18 


26827 


42 


s 





9 


1 1 


23 09.24, 


+42 55 38.3 


19.10 


17.39 


58448 


31 


s 


2 


10 


1 1 


23 10.32, 


+42 53 06.2 


18.49 


17.10 


33987 


41 


s 





1 1 


1 1 


23 11.96, 


+42 56 08.7 


20.40 


19.20 


58737 


78 


T 


*2 


12 


1 1 


23 14.37, 


+42 57 28.6 


21.00 


19.34 


90603 


72 


T 


*0 


13 


1 1 


23 15.00, 


+43 12 08.6 


19.24 


17.56 


24005 


30 


s 





14 


1 1 


23 15.76, 


+42 58 06.7 


23.01 


21.11 


56786 


150 


T 


2 


15 


1 1 


23 16.45, 


+42 53 50.3 


20.20 


18.05 


57516 


56 


T 


2 


16 


1 1 


23 18.74, 


+42 57 30.1 


20.38 


18.68 


73676 


88 


T 


*0 


17 


1 1 


23 19.68, 


+43 18 44.6 


21.96 


18.84 


133977 


49 


s 





18 


1 1 


23 20.92, 


+42 56 28.3 


21.40 


19.08 


57994 


100 


T 


2 


19 


1 1 


23 21.02, 


+42 57 12.5 


19.58 


17.91 


57551 


72 


T 


2 


20 


1 1 


23 21.69, 


+42 55 25.9 


20.22 


17.83 


57982 


52 


T 


2 


21 


11 


23 22.07, 


+42 54 58.6 


21.64 


19.22 


59445 


92 


T 


2 


22 


11 


23 22.07, 


+43 04 44.2 


20.81 


18.64 


59551 


83 


T 


1 


23 


11 


23 23.23, 


+42 54 27.7 


21.22 


19.03 


58431 


82 


T 


2 


24 


11 


23 23.70, 


+42 56 13.9 


20.29 


17.88 


58813 


38 


T 


2 


25 


11 


23 23.88, 


+42 58 25.7 


23.49 


19.78 


90451 


142 


T 


*0 


26 


11 


23 24.26, 


+43 09 38.1 


20.57 


18.12 


59406 


140 


T 


1 


27 


11 


23 24.64, 


+42 56 06.0 


19.46 


16.92 


57965 


38 


T + S 


2 


28 


11 


23 25.78, 


+42 51 34.7 


22.14 


20.45 


58742 


120 


T 


2 


29 


11 


23 25.80, 


+42 54 17.6 


19.04 


16.49 


57482 


52 


S 


2 


30 


11 


23 26.04, 


+43 05 10.5 


20.18 


17.86 


59031 


50 


T 


1 


31 


11 


23 27.35, 


+43 00 00.9 


20.28 


18.34 


35290 


70 


T 





32 


11 


23 27.80, 


+42 59 35.0 


18.51 


17.22 


17883 


26 


T + S 


*0 


33 


11 


23 28.33, 


+43 05 48.1 


21.23 


19.27 


57149 


174 


T 


1 


34 


11 


23 28.35, 


+42 53 49.1 


20.77 


18.52 


57489 


36 


T 


2 


35 


11 


23 28.43, 


+43 05 29.9 


20.49 


18.22 


57988 


61 


T 


1 



3. Analysis and Results 

3.1. Member selection 

To select cluster members out of 145 galaxies having redshifts, 
we follow a two steps procedure. First, we perform the ID 
adaptive-kernel method (hereafter DEDICA, Pisani 1993 and 
1996; see also Fadda et al. 1 19961 Girardi et al. U996b . We search 
for significant peaks in the velocity distribution at >99% c.l.. 
This procedure detects A1237+A1240 as a peak at z ~ 0.196 
populated by 95 galaxies considered as candidate cluster mem- 
bers (see Fig. [3j. Out of 50 non members, 24 and 26 are fore- 
ground and background galaxies, respectively. 

All the galaxies assigned to the cluster peak are analyzed in 
the second step which uses the combination of position and ve- 
locity information: the "shifting gapper" method by Fadda et al. 
(1996). This procedure rejects galaxies that are too far in ve- 
locity from the main body of galaxies within a fixed bin that 
shifts along the distance from the cluster center. The procedure 
is iterated until the number of cluster members converges to a 



stable value. Following Fadda et al. ( 1996) we use a gap of 1000 
km s _1 - in the cluster rest-frame - and a bin of 0.6 hzl Mpc, or 
large enough to include 15 galaxies. 

The choice of the center of A1240 is not obvious. No evi- 
dent dominant galaxy is present, rather there are some luminous 
galaxies. In particular, the two brightest ones, ID. 75 and ID. 56, 
lie in the southern and northern region of A 1240, respectively, 
and show comparable luminosity (hereafter BCGS and BCGN). 
The third one is located in the eastern region, but is > 0.5 R— 
magnitudes fainter than BCGS and BCGN (ID. Ill, hereafter 
BCGE). As for the cluster center, hereafter we assume the posi- 
tion of the centroid of the X-ray emission as recovered by David 
et al. JT9991 [R.A. = 1 l h 23 m 37 ?6, Dec.=+43°05'51" (J2000.0)] 
which lies between the two dominant galaxies. After the "shift- 
ing gapper" procedure we obtain a sample of 89 fiducial cluster 
members (see Fig.|4|i. 

We also check the result of alternative member selection pro- 
cedures. We apply the "shifting gapper" procedure adopting as 
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ID 



a, S (J2000) 



v 

(km 



) 



36 11 

37 11 

38 11 

39 11 

40 11 

41 11 

42 11 

43 11 

44 11 

45 11 

46 11 

47 11 

48 11 

49 11 

50 11 

51 11 

52 11 

53 11 

54 11 

55 11 

56 11 

57 11 

58 11 

59 11 

60 11 

61 11 

62 11 

63 11 

64 11 

65 11 

66 11 

67 11 

68 11 

69 11 

70 11 



23 28.81 
23 29.80 
23 29.97 
23 30.13 
23 30.17 
23 30.84 
23 31.07 
23 31.17 
23 31.20 
23 31.32 
23 32.10 
23 32.40 
23 32.68 
23 32.75 
23 33.08 
23 33.72 
23 33.84 
23 34.15 
23 34.67 
23 34.82 
23 35.26 
23 35.46 
23 35.48 
23 35.52 
23 35.60 
23 35.77 
23 36.00 
23 36.34 
23 36.56 
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cluster center the brightest cluster galaxy (BCGS). We find a 
cluster sample of 89 galaxies, 88 of which are in common with 
our above sample. In order to analyze the effect of a fully alterna- 
tive selection procedure we also apply the classical 3-cr clipping 
procedure by Yahil & Vidal (1977) on the whole sample of 145 
galaxies after a very rough cut in the velocity space, i.e. rejecting 
galaxies with velocities differing by more than 8000 km s _1 from 
the mean velocity. This classical procedure leads to a sample of 
90 galaxies, 89 of which forms our adopted sample. In conclu- 
sion, the sample of member galaxies we adopt in this work is 
quite robust against the member selection procedure. 

The galaxy distribution analyzed through the 2D DEDICA 
method clearly shows the presence of a southern external clump 
(see Fig. [5] see also Sect. 13 .7t . Gal et al. (2003 ) recovered a 
cluster in the same position from the digitized Second Palomar 
Observatory Sky Survey. We identify this galaxy clump with 
A 1237 which is likely to have a cluster redshift similar to that 
of A 1240 (cf. the magnitudes of their lOth-ranked cluster galax- 
ies, Abell et al. 1 19891 >. Notice that the center reported by Abell 
et al. is quite imprecise and lies on the southern border of the 
galaxy concentration we detect. 



We use the 2D DEDICA results, i.e. the peaks detected in 
the 2D galaxy distribution, to assign galaxies to different sub- 
clumps. The 2D DEDICA algorithm detects nine peaks, four 
of which are more significant than the 99% c.l.. The southern 
three peaks, only one of which is very significant, are assigned 
to A1237 (for a total of 27 members). The six northern peaks, 
three of which are very significant, are assigned to A 1240 (for a 
total of 62 members). This assignment is shown in Fig. [6]- left 
panel. 

As for A1240, the 2D DEDICA algorithm shows a clear bi- 
modal structure (see Fig. [5]) along the North-South direction. 
This bimodality is also shown in our analysis of photometric 
"likely" members in Sect. |3.7| and corresponds to the elongated 
hot gas distribution shown by the previous analyses of ROSAT 
data (David et al. 1999; Bonafede et al. 120091 ) and our analysis 
of Chandra data (see Sect.|4]). Therefore we decide to consider: a 
southern structure - hereafter A1240S - associated to the south- 
ern peak of A 1240 (the most significant in the whole DEDICA 
analysis); a northern structure - hereafter al240N - associated 
to the four northern peaks (two of which are very significant). 
In this way we assign 32 (27) members to A1240N (A1240S). 
A1240S and A1240N host the brightest and the second bright- 
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est galaxies, BCGS and BCGN, respectively. We consider sep- 
arately three galaxies belonging to a minor, eastern peak (here- 
after A1240E) since their assignation to A1240N or A1240S is 
not obvious. A1240E hosts the BCGE. The assignment of galax- 
ies within A 1240 is summarized as follow: we assign to A1240S 
the galaxies belonging to the southern peak (detected with a 
> 99.99% c.l.); to A1240N the galaxies belonging to the four, 
northern peaks (two of which are detected with a > 99.8% c.l.); 
to al240E the galaxies belonging to the eastern peak. 
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Fig. 3. Redshift galaxy distribution. The solid line histogram 
refers to the 95 galaxies assigned to the A1237+A1240 complex 
according to the DEDICA reconstruction method. 



O 



-2 - 



n 1 1 r 



h 1 1 r 



n 1 1 r 



do % ° 

o ® ° 

o o 

O x 



o o ° O 
<3>' 



O 



_l I l_ 



_L 



x 

_i I i i i_ 



12 3 

R [h" > Mpc] 

Fig. 4. Rest-frame velocity vs. projected clustercentric distance 
for the 95 galaxies in the main peak (Fig. 01 showing galax- 
ies detected as interlopers by our "shifting gapper" procedure 
(crosses). Circles indicate the 89 cluster members. 



3.2. The A1237+A1240 complex in the velocity space 

According to the ID Kolmogorov-Smirnov test (hereafter 
lDKS-test; see, e.g., Press et al 1992), there is no significant dif- 
ference between the velocity distributions of A 1237 and A 1240 
(see also Fig. [6]- right panel). This result suggests us to investi- 
gate the global velocity distribution of the complex. 



We analyze the velocity distribution of cluster galaxies (see 
Fig- 13 using a few tests where the null hypothesis is that the ve- 
locity distribution is a single Gaussian. We estimate three shape 
estimators: the kurtosis, the skewness, and the scaled tail index 
(see, e.g., Beers et al. |19911 l. According to the value of the skew- 
ness (-0.379) the velocity distribution is marginally asymmetric 
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Fig. 5. Spatial distribution on the sky and relative isodensity con- 
tour map of cluster members (A1237+A1240), obtained with the 
DEDICA method. The positions of the brightest galaxies are in- 
dicated (BCGN, BCGS, BCGE for A1240 and BCG1 and BCG2 
for A1237). The plot is centered on the cluster center defined in 
this paper as the X-ray center. The two relics are indicated in a 
schematic way. 
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Fig. 6. Spatial distribution on the sky of 89 cluster members. 
Left panel, open symbols indicate cluster members. In particular, 
squares, circles, rotated squares and triangles indicate A1240N, 
A1240S, A1240E and A1237. The crosses indicate emission line 
galaxies. Right panel: the larger the symbol, the smaller is the 
radial velocity. The six close northern circles indicate the high 
velocity galaxies we consider responsible of the apparent large 
global mean velocity of A1240N galaxies (see Sect. 13.41 and 
Fig.HOlin the following, too). 
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Fig. 7. Upper panel: velocity histogram for the whole cluster 
sample (thin line), A1240 (solid line) and A1237 (dashed line). 
The velocities of the three brightest galaxies of A 1240 are indi- 
cated (BCGN, BCGS and BCGE), as well as the two brightest 
galaxies of A1237 (BCG1 and BCG2). Lower panel: stripe den- 
sity plot where the arrows indicate the positions of the signifi- 
cant gaps. The main and minor weighted gaps are indicated by 
big and small arrows, respectively. 



differing from a Gaussian at the 90-95% c.l. (see Table 2 of Bird 
& Beers 1993). We also analyze the presence of "weighted gaps" 
in the velocity distribution. A weighted gap in the space of the 
ordered velocities is defined as the difference between two con- 
tiguous velocities, weighted by the location of these velocities 
with respect to the middle of the data (Wainer and Schacht 1978 ; 
Beers et al. |1991| l. We detect a strongly significant gap (at the 
> 99.9% c.l.) and five minor gaps (at the > 98% c.l.), see Fig. [7] 
- lower panel. The most important gap, very significant since 
it is located in the central region of the velocity distribution, 
separates the cluster into two subgroups of 37 and 52 galaxies. 
When comparing the 2D galaxy distributions of these subgroups 
we find no difference according to the 2D Kolmogorov-Smirnov 
test (Fasano & Franceschini[~ 



We also perform the 2D and 3D Kaye's mixture model 
(KMM) test (as implemented by Ashman et al. 1994) and com- 
pare the results to check the effect of the addition of the velocity 
information. The KMM algorithm fits a user-specified number 
of Gaussian distributions to a dataset and assesses the improve- 
ment of that fit over a single Gaussian and give an assignment of 
objects into groups. We use the A1237+A1240 galaxy assign- 
ment found by the 2D DEDICA analysis to determine the first 
guess when fitting two groups. The 2D KMM algorithm fits a 
two-group partition, at the > 99.9% c.l. according to the likeli- 
hood ratio test, leading to two groups of 66 and 23 galaxies. The 
addition of the velocity information in the KMM algorithm leads 
to the same group partition. 

Finally, we combine galaxy velocity and position informa- 
tion to compute the A-statistics devised by Dressier & Shectman 
( 1988 ). This test is sensitive to spatially compact subsystems that 
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have either an average velocity that differs from the cluster mean, 
or a velocity dispersion that differs from the global one, or both. 
We find no significant substructure. 

We conclude that, although the velocity distribution shows 
evidence of a complex structure, A 1237 and A 1240 are so sim- 
ilar in the velocity space that the velocity information is not 
useful to improve the galaxy assignment recovered from the 2D 
analysis (see Sect. 13. U . 

3.3. Global Kinematical properties 

As for the whole cluster complex, by applying the biweight es- 
timator to the 89 members (Beers et al. 1990), we compute a 
mean redshift of <z) = 0.1937+ 0.0003, i.e. <v> = (58273 ± 89) 
km s _1 . We estimate the LOS velocity dispersion, cry, by us- 
ing the biweight estimator and applying the cosmological cor- 
rection and the standard correction for velocity errors (Danese 
et al. 1 1980t . We obtain cry - 842+^ km s _1 , where errors are 
estimated through a bootstrap technique. 

The results obtained for the 62 members of A 1240 are: 
(z) = 0.1948+ 0.0004, i.e. (v) = (58408 + 111) km s" 1 , and 
<t v = 870^' km s . To evaluate the robustness of cry of A 1240 
we analyze the velocity dispersion profile (Fig. [8]). The integral 
profile is roughly fiat in the external cluster regions suggesting 
that the value of cry for A 1240 is quite robust. Figure [8] also 
shows that the (v) and cry profiles are not disturbed by the pres- 
ence of A1237. 

Table [2] lists the number of the member galaxies N g and the 
main kinematical properties of A 1237 and A 1240. 

Figure [9] compares the (v) and cry profiles of A1240 and 
A1237, where for A1240 we show separately the results for 
A1240N and A1240S. The value <v) of A 1237 is similar to 
that of A1240S, suggesting a continuity in the velocity field. 
However the value of cry of A 1237 is clearly lower than that of 
A1240S, suggesting that A1237 is really a less massive system. 



3.4. Internal structure of A1 240 

According to the lDKS-test there is no difference between the 
velocity distributions of A1240N and A1240S. The velocity dis- 
tribution of A 1240 shows signatures of non-Gaussianity similar 
to those of the whole cluster complex (e.g. a slight asymmetry 
and an important gap, see Pig - [7j . Like for the A1237 vs. A1240 
case, we find that A1240N and A1240S are so similar in the ve- 
locity space that the velocity information is not useful to improve 
the galaxy assignment. However, at the contrary of the A1237 vs. 
A 1240 case, A1240N and A1240S are spatially closer and likely 
strongly interacting (see Sect. [5j. This suggests that the galaxy 
assignment might be questionable and that we must devote more 
cure in determining the individual kinematical properties of the 
two subclumps. 

The value of global (v) of A1240N is nominally larger than 
that of A1240S, i.e. (v) = (58624 + 128) km s" 1 and (v> = 
(58091 + 195) km s _1 respectively, only at a ~ 2<x c.l. Moreover, 
looking at Fig. [9] the two cores of A1240N and A1240S seem to 
have similar (v): this is better shown in Fig[l0] where we directly 
compare the profiles of A1240N and A1240S. The large global 
(v) of A1240N is likely due to few high velocity galaxies in the 
extreme northern cluster regions (see the six galaxies shown as 
close circles in Fig. [6]- right panel; see also Fig. [TO}. Indeed, the 
recent merger of two subclumps may result in a plume, or arm, 
of outflying galaxies detected for their different velocity with re- 
spect to the cluster (see e.g. the case of Abell 3266; Quintana et 
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Fig. 8. Upper panel: rest-frame velocity (in units of 10 3 km s _1 ) 
vs. projected clustercentric distance for the 89 cluster galaxies. 
Circles and squares indicate members of A 1240 and A 1237. 
Central and lower panels: integral profiles of mean LOS velocity 
(in units of 10 3 km s~') and LOS velocity dispersion (in units of 
km s" 1 ) for the 62 members of A1240 (dots with the 68% c.l. er- 
rorbar) and for all the 89 cluster members (A1237+A1240, solid 
line). The mean and dispersion at a given (projected) radius from 
the cluster-center is estimated by considering all galaxies within 
that radius - the first value computed on the five galaxies closest 
to the center. In the lower panel, the horizontal lines represent the 
X-ray temperature with the respective 68% errors transformed 
in cry assuming the density-energy equipartition between gas 
and galaxies, i.e. yS spec = 1 where y6 spec = o% I \kT / pm p ) with 
p = 0.58 the mean molecular weight and m p the proton mass. In 
all the panels the center is the X-ray center of A 1240. 



al. 1 19991 and Flores et al. 120001 . Due to these difficulties in de- 
tecting a quantitative difference between A1240N and A1240S, 
we prefer to use the position and velocity of BCGN and BCGS as 
tracers of the two interveining clumps. In fact, dominant galax- 
ies trace the cluster substructures (Beers & Geller 1983) and are 
good tracers of interacting subclumps during a cluster merger, 
too (e.g., Boschin et al. 2006 ; Barrena et al. |2007at Boschin et 
al. l2M9l . 

The nominal value of global cry of A1240N is smaller than 
that of A1240S, cry = 709+g km s" 1 and cry = 991+^ 9 
km s _1 respectively. Using the brightest galaxies as centers the 
comparison of the respective cry profiles confirms this trend (see 
Fig[10]- right lower panel). Thus, although the nominal values 
of individual cry might be not fully reliable, we decide to adopt 
them for describing A1240N and A1240S. 

Table [2] summarizes the properties of A1240N and A1240S. 

3.5. Internal structure of A1 237 

As for A1237, it appears dominated by two bright galaxies, 
BCG1 and BCG2 (IDs. 29 and 27, see Table [B and the 2D 
DEDICA algorithm shows the presence of one significant peak 
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Table 2. Global properties of galaxy systems. 



System N g a, S < v > try R v „ M{< R^) 
(J2000) kins' 1 kms' 1 fo^'Mpc 10 14 M E 

A1237 27 112329.5 +425419 58021 ± 145 738+jg 1.6 6 + 2 

A1240 62 112337.6 + 430551 58408 + 111 870+'' 1.9-2.4 b 9 - 19 b 

A1240N" 32 112335.2 + 430832 - 58353 709™ 1.6 5 ±2 

A1240S a 27 112337.7 +430328 - 58817 991+if 2.2 14 ± 5 



a As center and mean velocity of this clump we list the position and velocity of the respective brightest galaxy. 

b The lower limit comes from the observed cr v . The upper limit is obtained when considering the bimodal structure of the cluster (see text). 
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Fig. 9. Upper panel: rest-frame velocity (in units of 10 3 km s~') 
vs. projected clustercentric distance for the members of A1240N 
(squares), A1240S (circles) and A1237 (crosses). In this and 
other panels the lower axis give the clustercentric distance of 
A1240N (A1240S) galaxies versus the negative (positive) axis 
from the A 1240 center, while the upper axis give the cluster- 
centric distance of A 1237 galaxies from the A 1237 center. The 
center of A 1240 is defined as the X-ray center. The center of 
A 1237 is defined as the position of the main peak in the 2D 
analysis. Central and lower panels: 68% error bands of integral 
profiles of LOS mean velocity (in units of 10 3 km s~') and LOS 
velocity dispersion (in units of km s~') for A1240N and A1240S 
(solid lines) and A1237 (dashed lines). The circles in the central 
panel indicate the positions and velocities of the two brightest 
galaxies of A1240 (BCGN and BCGS) and the two brightest 
galaxies of A1237 (BCG1 and BCG2). 
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Fig. 10. Upper panels: Error bands of integral profiles of LOS 
mean velocity for A1240N (dashed line) and A1240S (solid 
line). The distance is computed using as the center the position 
of the most important peak (on the left) and the position of the 
brightest cluster member (on the right). The horizontal lines indi- 
cate the velocities of BCGN and BCGS (dashed and solid lines, 
respectively). The close circles point out at which radius the six 
northern high velocity galaxies in A1240N are included to show 
their effect on the estimate of the mean velocity (see the text and 
Fig- [SI too). Lower panels: The same for integral profiles of the 
LOS velocity dispersion. 



tion shows evidence for a platykurtic behavior, but at a poorly 
significant level (at the - 95% c.L). 



(at the > 99% c.l.) and two minor peaks significant at the 97% 
and 89% c.l. (see the somewhat elongated structure of A1237 
in Fig. [5]). Both the analysis of SDSS and our photometric data 
show that A 1237 has only one well defined peak in the 2D 
galaxy distribution (see Sect. 13. 1\ . This last result is based on 
a much larger sample, thus we conclude that we have no evi- 
dence for a complex structure in A1237. The velocity distribu- 



3.6. Mass estimates 

Under the usual assumptions (cluster sphericity, dynamical equi- 
librium and galaxy distribution tracing the mass distribution), 
we can compute global virial quantities. Following the prescrip- 
tions of Girardi & Mezzetti (2001), the virial radius is /? v ; r = 
0.17 x o-y/H(z) ft™ Mpc (see their Eq. 1 after introducing the 
scaling with H(z); see also Eq. 8 of Carlberg et al. 119971 for 
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R200) and the virial mass (Limber & Mathews 1960; see also, 
e.g., Girardi et al. [1998) is: 

M = 3n/2 ■ o-yRpv/G - SPT. (1) 

The quantity SPT, the surface pressure term correction (The & 
White 119861 ), is assumed to be equal to the 20% of the mass 
since this is the typical value recovered for the mass computed 
within the virial radius in the literature when the data of many 
clusters are combined together to enlarge the statistics (Carlberg 
et al. 119971 Girardi et al. 11998b . The quantity R PV is a pro- 
jected radius (equal to two times the projected harmonic ra- 
dius). The value of Rpy depends on the size of the sampled re- 
gion and possibly on the quality of the spatial sampling (e.g., 
whether the cluster is uniformly sampled or not). It is also 
possible to use an alternative estimate of Rpy based on a pri- 
ori knowledge of the galaxy distribution (see the Appendix in 
Girardi et al. 1995). Following Girardi et al. (1998; see also 
Girardi & Mezzetti 120011) we can assume a King-like distri- 
bution with parameters typical of nearby/medium-redshift clus- 
ters: a core radius R c = 1/20 x R w \ r and a slope-parameter 
jSfit.gai = 0.8, i.e., the volume galaxy density at large radii 
scales as r~ 3iSr "-s al = r~ 2A . With these assumptions we can use 
the eq. A6 of Girardi et al. ( 119951 ) to estimate Rpy, see also 
eq. 13 of Girardi et al. 1998] for a useful approximation (i.e., 
R PV = 1.1897? vir [l+0.053(7? vir /#c)]/[l +0.117(7?™/^)]. Having 
assumed the galaxy distribution, the value of Rpy depends only 
on R wn . In this way our estimates of global virial quantities only 
depend on our estimate of cry. 

Through this procedure we obtain a mass estimate Mai237(< 
R wn = 1.63 fc^jMpc) = (6 ± 2) xlO 14 M and M A i 24 o(< 
R vk = 1.92 ft^Mpc) = (9 + 3) xlO 14 Kj\ M . 

When a cluster shows a strongly substructured appearance 
(e.g., a bimodal structure), the use of the global cry to com- 
pute the mass might be misleading (Girardi et al. 1 19971 and refs. 
therein). The true mass could be overestimated or underesti- 
mated depending on the angle of view of the cluster structure. 
When the two subclumps are aligned along the LOS, they can- 
not be clearly distinguished in their projection onto the sky, but 
they can appear as two peaks (less or more overlapped, depend- 
ing on their relative velocity) in the redshift distribution: in this 
case the mass estimated by the observed global cry is likely to be 
an overestimate of the true cluster mass (e.g., cf. Tabs. 7 and 8 of 
Girardi et al. 1998 ). When the two subclumps are aligned along 
an axis parallel to the plane of sky, they appear as two structures 
in their projection onto the sky but they cannot be distinguished 
in the redshift distribution (since their relative velocity has no 
component along the LOS direction). In this case the global ve- 
locity distribution is likely formed by the velocity distributions 
of the two clusters somewhat superimposed and the observed 
global cry does not take into account the existence of both the 
two subclumps. 

The A1240N+A1240S system is comparable to the second 
case discussed above: the two subclumps are distinguished in 
the sky, but lie at a similar z. Indeed our analysis of Sect. [5] 
will show that the axis of the A1240N+A1240S system is likely 
roughly parallel to the plane of sky. A more reliable mass esti- 
mate of A 1240 might be obtained adding the mass estimates of 
the two subclumps, each mass recovered from their respective 
cry (see Table©. We obtain a mass M A i240 ~ 1.9 xlO 15 M Q . 
Another possible approach is to consider the future, virialized 
A1240 cluster and its global properties. Since the cluster virial 
mass computed within the virial radius scales with ~ fx 3 ,, we 
expect that the cr v vir of the virialized A 1240 corresponds to 



(°"v,A1240N + °"V,A1240 S ) 1/3 ' Le ' °~V,vir ~ 1100 km S" 1 Well 

larger than the LOS cry measured on observed data. This cor- 
responds to a mass of Mai24o(< 7? v ; r = 2.4 /i^jMpc) ~ 1.9 
xlO 15 /i^q M e in good agreement with the above estimate. 

In conclusion, we estimate a mass range of Mai 240 = 
(0.9 - 1.9) xlO 15 h~\ M and of M a123 7 + ai240 = (1.5 - 2.5) 
xlO 15 h~l M Q . 

3.7. Analysis of photometric data 

The results of the 2D DEDICA method applied to the 89 clus- 
ter members are shown in Sect. 13. II However, our spectroscopic 
data suffer from magnitude incompleteness and the field around 
the cluster is not covered in an homogeneous way. To overcome 
these limits we use our photometric catalog. 

We select likely members on the basis of the B-R vs. R rela- 
tion, as already performed in some previous works of ours (e.g., 
Barrena et al. 2007b). To determine the relation we fix the slope 
according to Lopez-Cruz et al. (120041 see their Fig. 3) and apply 
the two-sigma-c lipping fitting procedure to the cluster members 
obtaining B-R = 3. 557 -0.069x7? for the red sequence of A1240 
galaxies. Figure [TTI shows B-R vs. R diagram for galaxies with 
available spectroscopy and the fitted line. 

Out of our photometric catalog we consider galaxies (objects 
with SExtractor stellar index < 0.9) lying within 0.25 mag of the 
relation. Following Visvanathan & Sandage (119771 ) the width of 
0.25 mag approximately corresponds to 2<x around the color- 
magnitude relation (see also Mazure et al. 1988 for a classical 
application to Coma cluster). The contour map for 370 likely 
cluster members having R < 20.5 shows the bimodal structure 
of A 1240 and the presence of A 1237, confirming the results of 
Sect. l3.1l (see Fig.fTST). Similar results are obtained with different 
magnitude cut-offs (e.g., R < 20 and R < 21). 

In this paper we also use public photometric data from the 
SDSS. In particular, we use r', /', z' magnitudes, already cor- 
rected for Galactic extinction and consider galaxies within a ra- 
dius of 30' from the cluster center. 

Following Boschin et al. (120081 see also Goto et al. 2002), 
out of SDSS photometric catalog we consider galaxies (here ob- 
jects with r' phototype parameter = 3) lying within +0.08 mag 
from the median values of r'-i' =0.47 and /'-z'=0.32 colors of 
the spectroscopically cluster members. The value of 0.08 mag 
is two times the typical scatter reported by Goto et al. ((2002) 
for the corresponding color-magnitude relations r'-i' vs. r' and 
i'—zf vs. r' . However, this member selection seems not enough 
conservative for the case of A 1240. In fact, using our spectro- 
scopic catalog, we notice that this color-color selection recog- 
nizes as "likely members" 21 out of the 56 non member galaxies 
(against the 1 1 out of 56 by using the color-magnitude selec- 
tion). Therefore we decide to use here a more conservative se- 
lection, i.e. considering only galaxies lying within +0.04 mag 
from the median values of r'-i' and i'-z' colors (see Fig. FPU with 
272 galaxies with r' < 20.8). The galaxy distribution shows the 
N-S elongation of A 1240 and the presence of A 1237. Moreover, 
another cluster is shown at ~ 25' South-East of A1240. This 
cluster was also detected by Gal et al. (2003 ) and Koester et al. 
(120071 ). who estimated a photometric z ~ 0.18-0.19. 

4. X-ray analysis of A1240 

The X-ray analysis of A 1240 is performed on the archival 
data of the Chandra ACIS-I observation 800407 (exposure ID 
#4961). The pointing has an exposure time of 52 ks. Data reduc- 



12 



R. Barrena et al.: Internal dynamics of Abell 1240 




R 

Fig. 11. B-R vs. R diagram for galaxies with available spec- 
troscopy is shown by large circles and squares (members and non 
members, respectively). Small points indicate galaxies found in 
our INT photometric sample, i.e. objects with a class star < 0.9. 
The solid line gives the best-fit color-magnitude relation as de- 
termined on spectroscopically confirmed member galaxies; the 
dashed lines are drawn at +0.25 mag from the color-magnitude 
relation. 

tion is performed using the package CIAC0 (Chandra Interactive 
Analysis of Observations) on chips 10, II, 12 and 13 (field of 

3 CIAO is freely available at http://asc.harvard.edu/ciao/ 
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Fig. 12. Spatial distribution on the sky and relative isodensity 
contour map of likely cluster members (applying the color- 
magnitude selection to the INT photometric catalog) with R < 
20.5, obtained with the DEDICA method. The positions of the 
brightest galaxies are indicated. The two relics are indicated in a 
schematic way. 
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Fig. 13. Spatial distribution on the sky and relative isodensity 
contour map of likely cluster members (applying the color-color 
selection to the SDSS photometric catalog) with r' < 20.8, ob- 
tained with the DEDICA method. The positions of the bright- 
est galaxies are indicated. The two relics are indicated in a 
schematic way. 

view ~ 17' x 17')- First, we remove events from the level 2 event 
list with a status not equal to zero and with grades one, five and 
seven. Then, we select all events with energy between 0.3 and 10 
keV. In addition, we clean bad offsets and examine the data, fil- 
tering out bad columns and removing times when the count rate 
exceeds three standard deviations from the mean count rate per 
3.3 s interval. We then clean the chips for flickering pixels, i.e., 
times where a pixel has events in two sequential 3.3 s intervals. 
The resulting exposure time for the reduced data is 51.3 ks. 

A quick look of the reduced image is sufficient to have 
an hint of the morphology of the extended X-ray emission of 
A1240. However, the cluster was centered exactly in the "cross" 
of the gaps among the four chips. To obtain a more precise re- 
sult, we have to correct the photon counts in the poorly exposed 
ACIS CCD gaps. First, we operate a binning of the reduced im- 
age. Then, we operate a soft smoothing and correct the photon 
counts with a division by an exposure map. The result is an im- 
age from which we extract the contour levels (soft photons in the 
energy range 0.5-2 keV) plotted in Fig.Q] 

The complex X-ray morphology of A 1240 is evident. In par- 
ticular, the central X-ray emission is clearly elongated in the N- 
S direction, i.e. the same direction defined by the two galaxy 
clumps A1240N and A1240S. Contour levels in Fig. Q] also re- 
veal a diffuse source ~6' SE of the cluster center. A visual in- 
spection of the INT /?-band image shows that this structure does 
not match with any evident galaxy concentration. We suspect it 
could be a distant background galaxy group. Unfortunately, a 
more quantitative morphological analysis is not trivial. In fact, 
the gaps among the ACIS chips make quite critical the compu- 
tation of the surface brightness distribution. Moreover, the dis- 
turbed morphology of the ICM does not justify the spherical 
symmetric assumption. This does not encourage the presentation 
of X-ray profiles. 

As for the spectral properties of the cluster X-ray photons, 
we compute a global estimate of the ICM temperature. The tem- 



R. Barrena et al.: Internal dynamics of Abell 1240 



13 




Fig. 14. 17'xl7' Chandra X-ray smoothed image (ID 4961) of 
A1240 in the energy band 0.5-2 keV (North at the top and East 
to the left). The central circle has a radius of ~2.58' (0.5 Mpc at 
the cluster redshift, see text). Black contours show the two radio 
relics. 



perature is computed from the X-ray spectrum of the cluster 
within a circular aperture of ~2.58' radius (0.5 Mpc at the clus- 
ter redshift; see Fig. [T4l around the center of the four ACIS 
chips. Fixing the absorbing galactic hydrogen column density 
at 1.99xl0 20 cm -2 , computed from the HI maps by Dickey & 
Lockman ( 1990), we fit a Raymond-Smith d 19771 1 spectrum us- 
ing the CIAO package Sherpa with a x 1 statistics and assuming 
a metal abundance of 0.3 in solar units. We find a best fitting 
temperature of Tx = 6.0 + 0.5 keV. 

We then search for temperature gradients by dividing the 0.5 
Mpc radius circle in four quadrants (chosen in order to avoid the 
gaps among the chips), as shown in Fig. [14] As expected for a 
merging cluster, there is some evidence that the ICM tempera- 
ture is not homogeneous, with the Western quadrant (7x=6.7^[2 
keV) somewhat hotter than the Northern (rx=5.1^ ( ' ) g keV), the 
Eastern (T x =5.2+J;^ keV) and Southern (r x =5.4+[^ keV) quad- 
rants. More detailed temperature and metallicity maps would be 
highly desirable to better describe the properties of the ICM, but 
the photon statistics is not good enough to this aim. In particu- 
lar, we do not provide any estimate of the ICM temperature in 
proximity of the two relics, where the X-ray surface brightness 
is too low to obtain any reliable measurement (see Fig.[l4l. 

5. Cluster structure: discussion 

Both optical and X-ray data indicate that A 1240 is a strongly 
substructured cluster, elongated in the N-S direction, the same 
of the axis of symmetry of the relics. We also detect two clumps, 
separated by 1 .2 Mpc, in the galaxy distribution. These ob- 
servational features suggest that we are looking at a cluster just 
forming through the merging of two main subclumps. The diffi- 
culty in separating the two subclumps in the velocity space and 
the small LOS velocity difference of the two BCGs suggest that 



the axis of the merger lies in the plane of the sky. The evidence 
of a very disturbed ICM distribution, somewhat displaced from 
the galaxy distribution, suggests that the merger is seen after the 
phase of the core passage, as shown by the results of the numer- 
ical simulation (e.g., Roettiger et al. 1 19971 their Figs. 7-14). 

As for the observational cluster structure, A 1240 is similar to 
Abell 3667 (see Roettiger et al. |1999l and refs. therein) where the 
optical and X-ray structures are elongated in a direction roughly 
similar to that of the axis of symmetry of the two radio relics 
which are separated by ~ 3 - 4 /i~ ( J Mpc. The two intervening 
galaxy subclumps are separated by ~ 1 \c\ Mpc with a small 
LOS velocity difference ~ 120 km s _1 between the dominant 
galaxies. Basic observational features of Abell 3667 were ex- 
plained with the "outgoing merger shocks" model, where shocks 
provide sites for diffusive shock acceleration of relativistic elec- 
trons causing the presence of the radio sources (Roettiger et al. 
[19991 . 

For A 1240 we have detected two galaxy subclumps in the 
N-S direction, the same direction of the elongation of the X- 
ray surface brightness and the axis of symmetry of the two ra- 
dio relics. The values of relevant parameters for the two-clumps 
system, as deduced from the BCGN and BCGS, are: the rela- 
tive LOS velocity in the rest-frame, V r f = 390 km s _1 , and the 
projected linear distance between the two clumps, D = 1.2 hji 
Mpc. The two roughly symmetric relics lie more externally, sep- 
arated by a projected linear distance ~ 2 h~ Q l Mpc (Bonafede et 
al. 120091 . 

We now use the above parameters and the mass of the system 
computed in a range M A i240 = (0.9 - 1.9) xlO 15 hz\ M Q to in- 
vestigate the relative dynamics of A1240N and A1240S. In par- 
ticular, we use different analytic approaches based on an energy 
integral formalism in the framework of locally flat spacetime and 
Newtonian gravity (e.g., Beers et al. |1982| l. 

First, we compute the Newtonian criterion for gravita- 
tional binding stated in terms of the observables as V^D < 
2GM sys sin 2 acosa, where a is the projection angle between the 
plane of the sky and the line connecting the centers of the two 
clumps. The thin curve in Fig. [15] separates the bound and un- 
bound regions according to the Newtonian criterion (above and 
below the curve, respectively). Considering the lower (upper) 
limit of Mai240, the N+S system is bound between 9° and 84° 
(6° and 89°); the corresponding probability, computed consider- 
ing the solid angles (i.e., jj cos a da), is 84% (89%). 

Then, we apply the analytical two-body model introduced 
by Beers et al. (1982) and Thompson ( 1982; see also Lubin et al. 
1998 for a more recent application). This model assumes radial 
orbits for the clumps with no shear or net rotation of the system. 
According to the boundary conditions usually considered, the 
clumps are assumed to start their evolution at time to — with 
separation do = 0, and now moving apart or coming together for 
the first time in their history. 

In the case of A1240N+A1240S system, where the first core 
passage has likely already occurred, we assume that the time 
to = with separation do — is the time of their core crossing 
and that we are looking at the cluster a time t after. To obtain 
an estimate of t we use the Mach number of the shock M ~ 3 
as recovered by Bonafede et al. (2009) from the radio spectral 
index. The Mach number is defined as M — v s /c s , where v s is 
the velocity of the shock and c s in the sound speed in the pre- 
shock gas (see e.g., Sarazin 2002 for a review). The value of c s , 
obtained from our estimate of Tx, i.e. c s ~ 995 km s _1 , leads 
to a value of v s ~ 3 x 10 3 km s _1 . Assuming the shock velocity 
as a constant, the shock covered a ~ 1 hzl Mpc scale (i.e. the 
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distance of the relics from the cluster center) in a time of ~ 0.3 
Gyrs. We assume this time as our estimate of t. Although the 
velocity of the shock is not constant, recent studies based on nu- 
merical simulations show how the variation of v s is much smaller 
than the variation of the relative velocity of the subclumps iden- 
tified with their dark matter components (see Fig. 4 of Springel 
& Farrar I20071 see Fig. 14 of Mastropietro & Burkert 2008), 
thus our rough estimate of t is acceptable at the first order of 
approximation. 

The bimodal model solution gives the total system mass A f sys 
as a function of a (e.g., Gregory & Thompson |1984| ). Figure 
compares the bimodal-model solutions with the observed mass 
of the system. The present solutions span the bound outgoing so- 
lutions (i.e. expanding), BO; the bound incoming solutions (i.e. 
collapsing), BI a and Bib; and the unbound outgoing solutions, 
UO. For the incoming case there are two solutions because of the 
ambiguity in the projection angle a. Both the BO and UO solu- 
tions are, in principle, consistent with the observed mass range. 
However, the BO solution is the more likely since the probability 
associated to the BO solution is much higher than that associated 
to the UO solution. In fact, we obtain that P(BO) and P(UO) are 
92% and 8%, respectively, where these probabilities are com- 
puted considering the solid angles (see above) and assuming that 
the region of M sys values is equally probable for individual so- 
lutions (see e.g. Barrena et al. l2007bV As for the projection an- 
gle we estimate a value of a ~ 10°. This small angle means 
that the real spatial distance between A 1240 subclumps is simi- 
lar to the projected one while the real, i.e. deprojected, velocity 
difference is ~ 2000 km s _1 , a quite reasonable value dur- 
ing cluster mergers (see, e.g. Springel & Farrar 2007 and refs. 
therein). Notice that the relative velocity between galaxy clumps 
is smaller than the shock velocity, i.e. the regime is not station- 
ary, but this is expected when comparing shock and collisionless 
components in numerical simulations (Springel & Farrar 2007 ; 
Mastropietro & Burkert [2008b . 

Since the merger of A 1240 clumps occurred largely in the 
plane of the sky, this likely explains the similarity between the 
observational features of A 1240 and Abell 3667. However, the 
case of A 1240 is made more complex by the presence of A 1237. 

We investigate the relative dynamics of A 1237 and A 1240 
with the same approach described above. The values of the rele- 
vant observable quantities for the two-clumps system are based 
on Table [2] i.e. LOS V rf = 320 km s 1 and D = 2.7 h~l 
Mpc. For the mass of the system we use the mass range com- 
puted in Sect. 13.61 Mai237+ai240 = (15 - 25) xlO 15 h~\ M . 
The Newtonian criterion for gravitational binding leads to a sys- 
tem bound between 8° and 89° (7° and 89°); the correspond- 
ing probability, computed considering the solid angles (i.e., 

f* 9 cos a da), is 87% (89%). 

In the case of A 1237 and A 1240 we have no evidence of a 
previous merger since the gas distribution of A 1240 (although 
very complex) does not show evidence of a peculiar displace- 
ment in direction of A1237. Therefore we assume that we are 
looking at A 1237 and A 1240 before their encounter. X-ray ob- 
servations of A 1237 would be useful to confirm this hypothe- 
sis. Under our assumption we can use the standard version of 
the analytical two-body model to study the A1237+A1240 sys- 
tem, where the clumps are assumed to start their evolution at 
time to — with separation do = 0, and are moving apart or 
coming together for the first time in their history. We are look- 
ing at the system at the time of t — 11.090 Gyrs, i.e. the age 
of the Universe at the system redshift. The solutions consistent 
with the observed mass span these cases: the bound and present 
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Fig. 15. System mass vs. projection angle for bound and un- 
bound solutions (thick solid and thick dashed curves, respec- 
tively) of the two-body model applied to A1240N and A1240S 
subsystems. Labels BI a and Bib indicate the bound and in- 
coming, i.e. collapsing solutions (thick solid curve). Label BO 
indicates the bound outgoing, i.e. expanding solutions (thick 
solid curve). Label UO indicates the unbound outgoing solu- 
tions (thick dashed curve). The horizontal lines give the range 
of observational values of the mass system. The thin dashed 
curve separates bound and unbound regions according to the 
Newtonian criterion (above and below the thin dashed curve, re- 
spectively). 

incoming solution (i.e. collapsing), BI a and Bib, and the bound- 
outgoing solution, BO (see Fig. [TBI). The associated probabilities 
are 61%, 31% and 8%, for BI a , BI b and BO respectively. The 
most probable solution is also the most interesting one since the 
BI a solution gives the same value of a ~ 10°already found for 
the A1240N+S merger. In this case, the direction of the clump 
velocities are consistent, too: at the present time A1240N has al- 
ready crossed A1240S, being now in front of A1240S and mov- 
ing outgoing; A1237 is lying behind A1240 and is infalling onto 
it along almost the same N-S direction. 

6. Conclusions 

Our results strongly support the conclusion of Bonafede et al. 
(2009), based on radio data, in favor of the "outgoing merger 
shocks" model for the double relics of A 1240. In fact, we detect 
the intervening merging subclumps and recover an acceptable 
model for the internal cluster dynamics. Moreover, our results 
also suggest that we are looking at the cluster accretion along a 
large scale structure filament. 

Our analysis shows how powerful is the study of the internal 
cluster dynamics through the analysis of kinematics of cluster 
galaxies. Other insights into A1240 might be recovered from a 
better knowledge of galaxy properties, e.g. spectral signatures 
of past activity which could be useful to determine the relevant 
time-scales (see e.g., Ferrari et al. 120031 Boschin et al. 2004) 
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and from deeper X-ray observations (see, e.g Vikhlinin et al. 
1200 II for the discover of a "cold front" in Abell 3667). Finally, 
we point out that A 1240, with its displacement between colli - 
sional (gas) and collisionless components (galaxies), is an excel- 
lent candidate to study the properties of the dark matter, in par- 
ticular its collisional cross section, as already performed in other 
merging clusters (Markevitch et al. 2004; Bradac et al. 2008). 
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